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Abstract. In a ongoing programme aimed at studying galaxy 
counterparts of Damped Lya Absorbers (DLAs) we have ob- 
tained high resolution deep I-band imaging data of the field 
around the z em = 0.457 BAL QSO Q21 12+059. In the liter- 
ature this QSO is listed to have a candidate DLA at z a \, s = 
0.2039 along the line of sight. After subtraction of the QSO 
Point Spread Function (PSF) we detect a galaxy centred on the 
position of Q21 12+059. To help answer whether this galaxy 
is the DLA or the QSO host galaxy we retrieved a GHRS 
spectrum of Q21 12+059 from the HST-archive. This spectrum 
shows that there is no Lya absorption line at z a & s = 0.2039. 
This fact in combination with the perfect alignment on the sky 
of the galaxy and Q21 12+059 lead us to the conclusion that the 
galaxy must be the host galaxy of Q21 12+059. 

The host galaxy of Q21 12+059 is bright (Mf s = -23.6), 
and has a radial profile well fitted by a modified Hubble + de 
Vaucouleurs profile with R c = 0.5kpc and R e = 3.6kpc. Our 
results are well in line with the conclusion of earlier work done 
at lower redshifts, that bright low redshift QSOs preferentially 
reside in luminous, elliptical galaxies. The host of Q21 12+059 
is however, despite it's brightness, very compact when com- 
pared to early type galaxies at lower redshifts. 

Key words: quasars - absorption lines, quasars - Q21 12+059, 
galaxies - active, galaxies - elliptical, lenticular, cD 



1. Introduction 

We have during the last several years been undertaking pro- 
grammes aimed at detecting optical counterparts of the most 
HI rich class of QSO absorption line systems, the Damped 



Send offprint requests to: J.U. Fynbo 

* Based on observations made with the Nordic Optical Telescope, 
operated on the island of La Palma jointly by Denmark, Finland, Ice- 
land, Norway, and Sweden, in the Spanish Observatorio del Roque de 
los Muchachos of the Instituto de Astrofisica de Canarias. 
** Based on observations with the NASA/ESA Hubble Space Tele- 
scope, obtained from the data Archive at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc. under NASA contract No. NAS5-26555. 
Correspondence to: jfynbo@eso.org 



Lya absorbers or DLAs for short (M0ller and Warren 1993, 
1998; Fynbo et al. 1999, 2000a,b). DLAs cause Lya absorp- 
tion lines with damping wings corresponding to HI column 
densities larger than 2 x 10 20 cto~ 2 . This very large column 
density absorption occurs in regions of self shielding, cool gas, 
i.e. where one would expect stars to form. Furthermore, they 
contain nearly all the neutral gas in the universe (Wolfe et al. 
1995). Hence DLAs at high redshift are prime candidates for 
being the progenitors of present day galaxies. 

At redshifts z < 1.65 DLAs are both rarer and observa- 
tionally more difficult to detect since Lya is not shifted above 
atmospheric cut-off and hence cannot be observed with ground 
based telescopes (Lanzetta et al. 1995; Jannuzi et al. 1998; Rao 
& Turnshek 2000). Le Brun et al. (1997) reported candidate 
galaxy counterparts of 7 such low redshift absorbers that are 
either confirmed DLAs or strong Mgll/Fell absorbers. They 
found that those candidate galaxy counterparts cover a wide 
range in morphological types. 

The z em = 0.457 BAL QSO Q2 112+059 (2000.0 coordi- 
nates RA 21 14 52.6, DEC +06 07 41) is radio quiet, and has 
optical magnitudes B = 15.52 and V = 15.77. Lanzetta et al. 
(1995) reported the detection of a candidate DLA at z a f, s = 
0.2039 in the IUE spectrum of Q21 12+059. In September 1996 
we performed an imaging study of Q21 12+059 in an attempt to 
detect the optical counterpart of this DLA. The imaging and 
data reduction is described in Sect. 2, in Sect. 3 we present 
spectroscopic data from the Hubble Data Archive and in Sect. 
4 and 5 we discuss our results. 

For easy comparison to earlier work on this subject, we 
shall here assume a Hubble constant of 50 km s _1 Mpc -1 and 

n=i.o. 



2. Imaging 

2.1. Observations and Data Reduction 

The observations were performed in 1996 with HiRAC (High 
Resolution Adaptive Camera) on the 2.56 m Nordic Optical 
Telescope under photometric and excellent seeing conditions. 
The CCD used was a lkxlk back-side illuminated thinned 
Tektronix with a pixel scale of 0.176 arcsec. The data were 
all obtained in the Bessel I filter. The journal of observations is 
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Fig. 1. Combined frame containing 973x1023 pixels 
(171 x 180 arcsec 2 ) surrounding the QSO Q21 12+059 (marked 
by an arrow). North is up and east to the left. Also marked are 
the field stars A-E which we use to check the stability of the 
PSF in Sect. 2.2. 

given in Table [l]. A total integration time of 3.25 hours was 
obtained. Also observed were 7 Landolt sequences (Landolt 
1992). 

The 39 individual frames were bias subtracted, flat-fielded 
and subsequently coadded using the optimal combination code 
described in M0ller and Warren (1993). 

In Fig. U we show the final combined image of 973 x 1023 
pixels (171x180 arcsec 2 ) surrounding Q21 12+059 which is 
marked with an arrow. Using the Landolt sequences we 
obtained the calibration from the instrumental system to 
Cousins I. The uncertainty in the zero-point is 0.02 mag. For 
Q21 12+059 we find I = 15.16±0.02, where the error is domi- 
nated by the zero-point error. The sky brightness was /.ij (sky) 
= 19.0 mag arcsec ~ 2 and the sky noise in the final combined 
image 26.8 mag arcsec ~ 2 . The full-width-at-half-maximum 
(fwhm) of point sources in the combined image is 0.69 arcsec. 



Table 1. Journal of NOT observations 
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Fig. 2. - a) The radial profiles of 5 bright field stars and of the 
QSO. The profiles were determined as the azimuthal average 
on a 9 by 9 sub-pixel grid. The profiles have been normalised 
to the same value in the centre (illustrated in the insert where 
we plot the central part of the fi\ profiles vs. log(R). The QSO 
clearly has excess extended flux compared to the stellar pro- 
files, b) Residuals after subtraction of the weighted mean stel- 
lar profile determined from the stars A-E. The residuals for the 
5 stars are to within the noise, but there is a clearly detected 
residual for Q21 12+059. 



Date 


Seeing FWHM 
(arcsec) 


Number of 
Frames 


Exp. Time 
(sec) 


1996 Sep 15 


0.6-0.8 


7 


2100 


1996 Sep 17 


0.5-0.7 


10 


3000 


1996 Sep 18 


0.6-0.8 


13 


3900 


1996 Sep20 


0.5-0.7 


4 


1200 


1996 Sep21 


0.8-0.9 


5 


1500 



There is no indication of variability of Q2 112+059 with 
amplitude above 0.01 mag on time-scales from tens of minutes 
to days. 
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2.2. PSF quality and stability test 

In order to search for faint sources near the line of sight to 
the bright Q21 12+059 a good PSF-subtraction is mandatory. To 
test the stability of the PSF we selected 5 bright stars distributed 
over the full field. These stars are marked A-E in Fig. |l| We 
then determined the radial profiles of the field stars by aver- 
aging over the azimuth angle on a 9 by 9 sub-pixel grid. In 
Fig. |2|a we have plotted the radial profiles determined in this 
way. Also plotted is the radial profile of Q21 12+059. All pro- 
files were normalised to the same level in the centre (shown in 
detail in the inserted logarithmic plot in Fig. ||a). The surface 
brightness (pi) zero-point in Fig. |2p was chosen to be that of 
the Q21 12+059 profile. It is seen that, i) the 5 bright stars have 
radial profiles (full drawn lines) that are identical to within the 
thickness of the lines out to about 2.8 arcsec, and ii) the ra- 
dial profile of Q21 12+059 (short dash) shows excess flux over 
that of the field stars. To quantify the excess we plot in Fig. ||b 
the residual surface brightness after subtraction of the weighted 
mean stellar profile determined from the field stars A - E. At 
radial distances ranging from 1.4 to 3.6 arcsec from the cen- 
troid of the QSO, the excess surface brightness is 0.08-0.14 
mag. In the same range the scatter of the stellar PSFs is about a 
factor 10-20 smaller. At radial distances larger than 3.6 arcsec 
the profiles are dominated by sky noise. 

We end this section concluding that the PSF is very stable 
across the field and that the profile of Q2 112+059 has signifi- 
cant excess flux over that of a PSF. 

2.3. 2D PSF-subtraction 

To examine the nature of the excess flux we proceeded with a 
two dimensional PSF subtraction using DAOPHOT II (Stetson 
1987, 1999). For the determination of the 2D PSF we chose to 
use star A, which is the closest to Q2 112+059, but as demon- 
strated above our results are not dependent on this specific 
choice of PSF star. When we subtracted the PSF of the QSO us- 
ing the scaling of the PSF automatically chosen by DAOPHOT 
II this produced a residual consisting of an extended positive 
residual surrounding a negative "hole" in which the PSF had 
been over-subtracted. This is the typical residual one will see 
when a point source and an extended object are superimposed 
on roughly the same position (see e.g. Fig. 2b where the PSF 
subtracted QSO profile is in the centre but 0.13 mag 3 arc- 
sec out). To obtain a image of the underlying extended object 
without the "hole" from the PSF subtraction, we must impose 
some additional condition. A reasonable assumption is that the 
extended object is "flat" across the centroid of the QSO, so we 
shall use this as the extra constraint to break the degeneracy. 

Fig. U shows the result of this final two dimensional PSF- 
subtraction. In the upper left panel is shown 100x100 pixels 
( 17.6 x 17.6 arcsec 2 ) around Q21 12+059. The upper right panel 
shows the same area after the final subtraction of the PSF. The 
residual image in the upper right frame shows a fairly large ex- 
tended object (marked S), presumably a galaxy. The centroid of 
this object coincides with that of the QSO to within 0.4 pixel 
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Fig.3. Upper left : 100x100 pixels (17.6x17.6 arcsec 2 ) sur- 
rounding Q21 12+059. Upper right : Same as left frame, but 
after subtraction of the 2D PSF. Lower panels : For compari- 
son we show the field star marked C in Fig. 1 before and after 
PSF-subtraction. The images have all been smoothed by a 3 x 3 
box car filter. The contour levels are -4, 4, 6, 8, 12, 20, 38,... in 
terms of sigma of the sky in the unfiltered image, with the neg- 
ative contour being dotted. 

(0.07 arcsec). The isophotes of the residual show no clear ev- 
idence for ellipticity, rather they are consistent with e = 0. 
There are two additional objects that are both consistent with 
being point sources 5 and 7 arcsec east of S, and patches of 
faint extended fuzz both east and west of the object S. In the 
lower panels of Fig. |] we show for comparison the field star 
C (0.35 mag fainter than Q21 12+059) before and after PSF- 
subtraction. 

3. Spectroscopy from the HDA 

Our original aim with the observations was to detect the op- 
tical counterpart of the DLA, but the perfect alignment of the 
galaxy S with the QSO reported above suggests that we prob- 
ably have detected the QSO host galaxy rather than the galaxy 
causing the reported candidate DLA absorption. To help clar- 
ify this we searched the Hubble Data Archive (HDA) for spec- 
troscopic data of Q21 12+059. In this section we describe the 
results of this search. 

3.1. Retrieval of the spectra and Data Reductions 

We found two data sets suitable for this study. Goddard High 
Resolution Spectrograph (GHRS) spectra over the region from 
1330 A to 1610 A (grism G140L) was obtained in 1995 (PI 
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Fig. 4. The 80A region around Lya at z=0.2039 of the normal- 
ized and binned spectrum of Q21 12+059. Also shown is the 
la error per 0.72A bin (lower curve). The dotted line shows 
a damped Lya absorption line with a H I column density of 
2x 10 20 cm~ 2 at a redshift of z=0.2039. The presence of such a 
line a excluded at the 29<r level. 



Lanzetta) and in 1992 spectra of the region from 1610 A to 
3250 A was obtained with the Faint Object Spectrograph (FOS) 
(PI Bahcall). Only the FOS-spectra has been published (Jan- 
nuzietal. 1998). 

We retrieved the calibrated GHRS spectrum from the 
archive. Hence, the data reduction was limited to extracting 
the 1 dimensional spectrum and 1 a error spectra from the cal- 
ibrated science data files and using the calibrated wavelength 
solution files. 




(R in arcsec) 



3.2. Results 

Fig. H shows the spectrum of Q21 12+059 covering the region 
from 1425 A to 1505 A. Also shown is the la error per 0.72 A 
bin (lower curve). 

The presence of a damped Lya absorption line in the Lya 
forest of Q21 12+059 at z abs = 0.2039 is excluded at the 29cx 
confidence level. The FOS spectrum shows absorption due to 
gas in the Milky Way and intervening metal line systems at 
redshifts z=0.370 and z=0.418 as well as metal line systems 
associated with Q21 12+059, but no absorption lines at the red- 
shift z=0.2039 (Jannuzi private communication; Turnshek et al. 
in prep.). 

As discussed by Rao & Turnshek (2000) only 2 of 14 can- 
didate DLA systems from the IUE survey of Lanzetta et al. 
(1995) have been confirmed. Q21 12+059 adds one to the list of 
non-confirmed candidates. 



Fig. 5. Upper panel: The crosses show the surface brightness 
for the host galaxy determined by azimuthal averaging on a 
9 by 9 sub-pixel grid plotted against R 1 / 4 (in arcsec). Also 
plotted (full drawn lines) are 6 exponential profiles with scale 
lengths 0. 1-0.6 arcsec convolved with the PSF and azimuthally 
averaged in the same way as the data. As seen no exponential 
model provides a good fit to the data. Lower panel.Same as the 
upper panel except here we have overplotted the radial profile 
of a model consisting of an outer R 1 / 4 profile and a modified 
Hubble profile in the core (see text for details). As seen this 
model match the data well. 

4. Nature of the host galaxy 

As there clearly is no DLA absorber in the sightline towards 
Q21 12+059, and as the centroid of the extended object we have 
detected is closely aligned with the position of the QSO (to 
within 0.07 arcsec), we conclude that we have detected the host 
galaxy of Q21 12+059. We now return to a discussion of what 
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can be learned about the nature of this galaxy from our deep 
imaging. 

An attempt to fit a 2D galaxy model to the 2D PSF sub- 
tracted image shown in Fig 3 was unsuccessful as the signal- 
to-noise of the image was too poor to provide meaningful con- 
straints. However, the alignment of the QSO with the centroid 
of its host, and the fact that there is no evidence in the 2D 
residual image (Fig. g) for any elongation of the extended ob- 
ject, suggest that we could significantly lower the number of 
fitting parameters, and at the same time increase the signal- 
to-noise of the data to fit, if we returned to the ID regime. 
Again, as in the 2D case, we are left with the degenerate prob- 
lem of deciding how to scale the PSF we subtract from the data. 
Galaxy cores are known often to have central spikes (see e.g. 
M0ller et al. 1995) which could be due either to central point 
sources (AGN), to central density cusps or (which is the most 
likely) to a combination of both (for a discussion see Stiavelli 
et al. 1993). In the case of Q21 12+059 there clearly is a central 
AGN (the QSO), which is the point source we wish to separate 
out. Therefore, as in M0ller et al. (1995) we break the degener- 
acy by imposing the condition that the central profile must be 
isothermal in the sense that the galaxy core is fitted by a mod- 
ified Hubble profile. This condition closely corresponds to the 
requirement we used for the 2D PSF subtraction, namely that 
the residual must be "flat" or "smooth" in the centre. With this 
assumption we effectively assign the maximum possible flux to 
the central point-source, and therefore obtain a lower limit to 
the flux of stellar origin. However, the central density cusps of 
local early-type galaxies do not contain any significant fraction 
of the total stellar light (1-10% of the flux measured inside the 
galaxy core, M0ller et al. 1995), so the resulting uncertainty on 
the flux from the total stellar component is insignificant. 

For the ID PSF we again used the weighted mean of the 5 
field stars (all field stars were given weights corresponding to 
their flux, except star E which was given a lower weight be- 
cause of the nearby saturated star which caused variable sky 
across the profile). In Fig. ||a and Fig. |]b we plot the host 
surface brightness per square arcsec (/i/) against R 1 / 4 . Super- 
imposed on the data are profiles corresponding to a range of 
galaxy models. In Fig. ||a we have superimposed six exponen- 
tial disc profiles 

I{R) = J x exp(-R/R s ) 

with scale lengths R s (0.1-0.6 arcsec), and in Fig. |5|b the radial 
profile of a model with an R 1 / 4 outer profile 

I(R) = J x exp(-7.67 x {{R/R e f- 2b - 1)), 

and a modified Hubble profile 

I(R) oc {l + (R/R c ) 2 )-\ 

in the core. The models shown were produced on a 2D sub- 
pixel grid and convolved with the observed 2D PSF. The final 
ID model profiles were then determined as azimuthal averages 
using the same procedure as for the data. 

It is seen that the radial profile of the host is noisy in the 
central 0.5 arcsec where the errors from the PSF subtraction are 



Table 2. Size, absolute magnitude and redshift of the host of 
Q21 12+059, and of Radio Galaxies (RG), Radio Load Quasars 
(RLQ) and Radio Quiet Quasars (RQQ) in the work by McLure 
etal. (1999). 

object R e ,rms ^gal^QSO M ga i , rms z 
kpc 

Q2 112+059 16 51)35 =231 0457 

RG 14.9 , 5.4 25-3000 - 23.65 , 0.50 0.174-0.244 

RLQ 11.2,4.5 0.17-1.3 -23.89,0.29 0.173-0.258 

RQQ 9.8,4.8 0.06-5 -23.27,0.52 0.115-0.239 



large and at radii larger than 3.5 arcsec where the uncertainty 
in the measurement of the sky-level becomes dominant (sky- 
level errors are shown as short dash curves in Fig. ||). At radii 
between 0.5 and 3.5 arcsec the profile is very well determined. 

Fig. [|a shows that it is not possible to obtain an acceptable 
fit with an exponential disc profile for any scale length. In con- 
trast, as seen in Fig. |]b, a model with an R 1 / 4 outer profile with 
an effective radius of R e — 0.53 arcsec (3.6 kpc) and a modi- 
fied Hubble profile with a core radius of R c = 0.07 arcsec (0.5 
kpc) in the core provides an excellent fit. 

5. Summary and Discussion 

Based on an IUE spectrum of Q21 12+059 Lanzetta et al. 
(1995) reported a candidate damped Lya line at z=0.2039. Sub- 
sequent HST spectroscopy (Fig. has not confirmed this. This 
result is important in its own right as DLAs are very rare at 
z < 1. In the sample of DLAs presented in Wolfe et al. (1995) 
there are 10 DLAs at z<1.0. Of the 10 only 4 are at z<0.5, and 
one of those was the candidate DLA towards Q21 12+059. 

Having ruled out the presence of an intervening DLA at 
z = 0.2039, we find that we are forced to interpret the ellip- 
tical galaxy we detect under the PSF of Q2 112+059 as the 
host galaxy of the QSO. The host galaxy has an I-band mag- 
nitude of 18.8 (total magnitude determined from the model 
fit), and an effective radius of 3.6 kpc. 1=18.8 corresponds to 
MJ hs = —23.6 (no k-correction) in the assumed cosmology 
(Sect. 1). 

McLure et al. (1999) find that both radio galaxies (RGs), 
radio loud QSOs (RLQs) and bright radio quiet QSOs (RQQs) 
preferentially are hosted by massive, bright elliptical galaxies. 
Our detection of a bright elliptical galaxy hosting the radio 
quiet Q2 112+059 is consistent with this result. In Table |]we 
compare our results to those of McLure et al. (1999) and find 
that the host galaxy of Q2 112+059 is in the bright end of the 
absolute magnitude distribution, but has a smaller effective ra- 
dius than any RQQ host galaxy in their sample. 

Schneider et al. (1983) studied the brightest cluster galax- 
ies in 83 Abell clusters (with redshifts in the same range as 
the McLure et al. samples). They found a tight correlation be- 
tween the core radii and effective radii. The host galaxy of 
Q21 12+059 falls nicely on this correlation, but is smaller than 
any of the brightest cluster galaxies of the Schneider et al. sam- 
ple. 
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Compared to the work on local early type galaxies by 
Kormendy (1977), the effective radius of the host galaxy of 
Q21 12+059 falls in the overlapping region between "compact" 
and "normal" early type galaxies. 

Finally, since the host galaxy we have detected is a spheroid 
we can estimate the mass of the black hole using the correla- 
tions derived by Magorrian et al. (1998) assuming that the cor- 
relation found at z=0 is also valid at z=0.457. We derive a black 
hole mass of about 9x 10 9 Mq. 
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